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摘要 : 应 用 随机 扩 增 多 态 性 DNA (RAPD) 技术 及 细胞 色素 bp (Cyt 2 ) 基因 序列 ， 分 析 了 长 江口 及 其 南部 邻 
近 地 区 大 弹 涂 鱼 (Boleophthalmus pecinirostris) 群体 遗传 结构 及 种 群 历史 ， 结 果 显 示 这 一 区 域 大 弹 涂 鱼 群 体 的 中 
传 多 样 性 水 平 很 高 。(1) RAPD 分 析 结 果 : 从 100 个 10 碱 基 随机 引物 中 筛选 出 30 个 多 态 性 引物 ， 对 长 江口 九段 
沙 湿地 、 浙 江 舟 山 定 海 、 福 建 霞 浦 群体 各 20 个 体 进行 RAPD 分 析 ，3 个 群体 分 别 获 得 236、270、274 条 带 ， 共 301 
个 有 效 位 点 。3 群体 多 态 位 点 比例 (P)、Nei 基 因 多 样 性 指数 (万 )、Shannon 多 样 性 指数 (7 分 别 在 93.02% 一 96.35%、 
0.3890 一 0.4219、0.5618 一 -0.6044。 分 子 方差 分 析 (AMOVA) 表明 ， 大 弹 涂 鱼 群 体 间 存在 显著 的 遗传 分 化 (FF.= 
0.03004-0.03254，P<0.05)， 但 仅 3.2% 的 遗传 变异 来 自 群 体 间 。(2〉 Cyt 2 基因 序列 分 析 结 果 : 从 采 目 九段 沙 、 定 
海 、 暇 浦 及 浙江 荡 溪 的 42 尾 大 弹 涂 鱼 样本 中 ， 共 获得 33 个 Cyt 2 基因 单 倍 型 〈 序 列 长 1141bp)。4 群体 的 平均 单 
倍 型 多 样 性 〈71)、 核 甘酸 多 样 性 (x) 分 别 为 0.9814、0.0048。4 群体 间 的 遗传 分 化 指数 下 ,为 0.00043 一 0.07814， 
仅 1.93% 的 变异 来 目 群 体 间 (AMOVA 分 析 )， 而 基因 交流 值 却 达 14.50 一 30.79， 群 体 间 K2-P 遗 传 距离 为 0.0040 一 
0.0056， 从 而 显示 大 弹 涂 鱼 群 体 间 没有 发 生 明 显 的 地 理 分 化 。 以 Cyt p 基 因 序 列 构建 的 NJ 树 揭示 4 个 群体 的 个 体 组 
成 2 个 谱系 ， 但 这 2 个 谱系 与 地 理 分 布 并 不 相关 。 中 性 检验 、 错 配 分 析 和 网 络 亲 缘 关 系 分 析 皆 表明 大 弹 涂 鱼 群 体 
有 过 种 群 扩张 ， 扩 张 时 间 约 在 更 新 世 末 期 的 0.057 一 0.023 百 万 年 前 。 
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Abstract: In order to analysis the genetic structure and demographic history of Boleophthalmus pectinirostris, the 
random amplified polymorphic DNA (RAPD) technology and mitochondrial cytochrome b (Cyt b) gene sequences were 
used in the study. The results showed that there was high level genetic diversity lin the populations of B. pecinirostris In 
Yangtze River estuary and 1ts southern wet-land adjacent regions. (1) RAPD was conducted to analysis the genetic 
variation of B. pecinirostris in Jiuduansha (JDS), Dinghai (DH) and Xiapu (XP ) for 20 individuals per population, with 
30 polymorphic 10-based random primers selected from 100 ones. In total, 236, 270 and 274 bands were obtained for 
JDS, DH and XP, respectively and there was 301useful loci. The range of proportion of polymorphic locl (P), Nel’s gene 
diversity (HFH) and Shannon diversity index ( ) was: 93.02%—96.3%, 0.3890-0.4219, 0.5618-—0.6044, respectively. The 
results of AMOVA showed significantly genetic differentiation (F,=0.03004-0.03254, P<0.05) but there was only 3.2% 
genetic variation among populations. (2) Results of Cyt b gene: the total Cyt b gene length was 1 141bp and 33 
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haplotypes were obtained from 42 indlividuals collecting in 4 different regions, which included JDS, DH, XP and CIXl 
(CX). The average haplotype diversity (h) and nucleotide diversity (x) was 0.9184, 0.0048, respectively. The high level of 
Nn values (14.50—30.79) and low Fs values (0.0043-0.07814), the average K2-P distances (0.0040—0.0056), as well as 
the AMOVA analysis (total F,=0.0193, P>0.05) all revealed that the B. pecinirostris had no obviously geographical 
differentiation. However, the phylogenetic tree estimated by nelghbor-jolning (NJ) method indicated that the ndividuals 


of B. pecinirostris formed two monophyletic lineages, and the two lineages were not correlated with geographical 
populations. Neutrality tests, mismatch distribution and network analysis suggested a late Pleistocene population 


expansion for B. pecinirostris (0.037-0.023 million years ago). 
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大 弹 涂 外 (Boleophthalmus pectinirostris) ,俗称 
化 跳 、 星 跳 、 跳 鱼 等 ， 隶属 于 饰 形 目 (Perciformes) 
贫 虎 角 科 (Gobiidae)， 为 上 暧 水 广 温 广 盐 性 海洋 色 类 ， 
生活 在 河口 及 漳 间 市 的 淤泥 滩 深 及 红 树 林 区 ， 在 我 
国 主要 分 布 于 江苏 、 上 海 、 浙 江 、 福 建 、 广 东 、 广 
西 和 台湾 沿海 。 大 弹 涂 鱼 不 仪 是 深 受 环保 专家 关注 
的 滩涂 湿地 指示 生物 ， 而 且 是 我 国 南 方 最 其 利用 价 
值 的 滩涂 人 鱼 关 之 一 。 日 本 和 我 国 台 湾 目 上 世纪 60 
年 代 就 开始 养殖 , 我 国 东 南 沿 海 的 养殖 则 兴起 于 20 
世纪 90 年 代 ， 已 成 为 当地 名 特 优 养殖 品种 。 

近年 来 ， 由 于 大 规模 的 海 导线 开 友 、 围 海造 地 
以 及 海 症 污染 ， 大 幅度 地 减少 并 割裂 了 大 弹 浴 鱼 奖 
的 栖 明 地 ， 过 上 度 的 捕捞 也 造成 了 其 野生 种 群 数 量 的 
严重 下 降 。 研 发 规模 化 育苗 技术 并 有 效 养 护 野生 种 
群 ， 已 成 为 大 弹 涂 鱼 产 业 化 开发 的 重要 环节 ， 而 欣 
明 其 种 群 跟 传 结构 是 实施 这 一 环节 的 前 提 。 目 前， 
有 关 大 弹 涂 鱼 的 形态 特征 、 行 为 习性 、 养 殖 繁育 技 
术 、 性 腺 发 育 等 已 有 许多 报道 (Zhang et al, 2006 )， 
但 在 种 群 组 成 和 遗传 结构 等 方面 的 研究 却 很 少 (Jin 
et al, 2004; 上 Kanemorl et al, 2006; Zhang et al, 2000)， 
因而 开展 此 方面 的 研究 很 有 必要 。 

随机 扩 增 多 态 性 DNA 〈RAPD) 是 动物 遗传 
多 样 性 研究 的 单 用 技术 ， 有 具有 安全 、 人 快捷 、 人 简便 、 
高 效 等 特点 (Liu et al, 2004; Matoso et al 2004)。 男 
外 一 类 分 子 标 记 ， 线 粒 体 DNA CmtDNA) 以 严格 
的 母系 遗传 和 较 快 的 进化 速率 ， 已 成 为 一 种 探讨 物 
种 系统 发 生 、 种 内 韶 传 分 化 的 有 效 遗 传 标记 
(Verheyen et al, 2003; Wang et al, 2006; Zhu et al, 
2007)， 而 细胞 色素 (Cyt 5) 是 mtDNA 中 应 用 较 
多 的 分 析 鱼 类 种 群 遗 传 结构 的 分 子 片 段 (Han et al 
2007; Santos et al, 2006; Wang et al, 2006)。 长 江口 及 
其 邻近 的 沿 尾 地 区 是 我 国人 类 活动 干预 最 强烈 的 
地 区 ， 建 立 有 国家 级 梨 明 东 滩 和 九段 沙 湿 地 目 然 保 
护 区 ， 大 弹 深 旬 是 这 些 保护 区 最 具 标 志 意 义 的 滩涂 










































































鱼 类 之 一 。 本 研究 旨 在 通过 RAPD 技术 分 析 与 Cyt。 
基因 全 序列 的 测定 ， 了 人 解 这 些 地 区 大 弹 涂 鱼 天 然 群 
体 的 遗传 结构 与 发 生 历 史 ， 为 其 资源 保育 与 合理 利 
用 提供 基础 资料 。 


1 材料 与 方法 


1.1 样本 采集 与 选用 
通过 当地 渔民 议 章 网 具 抓 捕 ， 于 2005 年 5 月 

一 2006 年 8 月 采集 到 一 批 不 同 地 方 的 大 弹 涂 鱼 天 然 
群体 样本 。 其 中 ， 任 取 九 段 沙 湿 地 目 然 保护 区 《从 
称 JDS)、 浙 江 舟 山 定海 (人 简称 DH)、 福 建 霞 浦 县 
官 井 洋 (人 简 称 XP) 样本 各 20 尾 , 用 于 RAPD 分 析 ; 
任 取 JDS 样本 11 尾 、XP 样本 10 尾 、DH 样本 10 
尾 、 浙 江 英 溪 样 本 11 尾 ( 简 称 CX)， 用 于 Cyt 5b 基 
因 序 列 测定 。 另 选取 采 目 舟山 定海 的 弹 浴 全 

( Periophthalmus modestus ) 、 大 鱼 弹 涂 鱼 〈 己 
magnuspinnatus) (Wang et al 2006) 及 福建 霞浦 的 青 
弹 涂 鱼 (Scartelqos viridis) 各 1 尾 ， 测 定 其 Cyt 2 
基因 序列 作为 构建 系统 发 育 树 的 外 类 和 群 。 所 有 样本 
用 95% 乙 酬 回 定 ， 一 周 后 保存 于 70% 乙 本 ， 首 4C 
冰箱 备用 。 取 过 样 的 标本 保存 于 上 海 海 洋 大 学 旬 类 
标本 室 。 
1.2 基因 组 DNA 提取 

基因 组 DNA 提 取 采 用 第 规 的 “ 酚 - 氯 仿 ” 法 ， 参 

照 Liu et al (2004) 的 方法 并 做 适当 修改 。 具体 操作 如 
下 : 取 0.1g 左 右 背 鳍 或 尾鳍 ， 切 人 雄 ， 加 入 49S$uL 裂 
解 绥 冲 液 (Sambrook et al, 2002), 54L 生 日 酶 K( 浓 
度 为 20 mg/mL)，55 一 56'C 恒 温水 浴 锅 消化 2 一 3h 

(也 可 放置 过 夜 , 约 12h)。 加 入 等 体积 饱和 酚 ， 绥 
慢 转 动 20min，12000r 离 心 10min。 吸 取 上 清 液 ， 
再 加 入 等 体积 的 氯仿 - 异 戊 醇 混 合 液 (24 : 1)， 绥 
慢 转 动 Smin。 取 上 清 液 , 加 入 2 倍 体积 的 预 冷 无 水 
乙醇 沉淀 DNA (-20C)， 用 70% 乙 醇 洗 涤 1 次 , 干 
燥 ， 加 入 300 一 5004L 的 ddH20 及 1 一 2nL RNA 酶 ， 
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4C 冰 箱 保存 备用 。 待 DNA 深 解 后 ， 取 5 4L 深 液 ， 
1% 琼 脂 糖 凝 胶 电泳 检 伍 DNA 提 取 和 情况 。 

1.3 RAPD-PCR 扩 增 

10 碱 基 随机 引物 100 个 (引物 编写 为 S1-S60， 

S101-S140)， 购 自 上 海 生 工 生物 工程 技术 服务 有 限 
公司 。 参照 说 明 书 配置 浓 虐 5pmol/L, 一 20"C 冰 箱 保 
存 备用 。 扩 增 条 件 参照 Liu et al (2004) 的 方法 。 取 
10mL 扩 增 产物 经 1% 进 口 琢 脂 糖 凝 胶 〈 含 EB 0.5 
mg/mL) 电泳 2h 后 ,用 Gene-Genius 凝 胶 成 像 系 统 
拍照 、 记 录 和 分 析 。 

1.4 Cytb 基因 的 PCR 扩 增 、 纯 化 及 测序 

扩 增 Cyt b 基因 使 用 的 特 寞 性 引物 为 : L14724 

(5’ -GAC TTG AAA AAC CAC CGT TG-3' ) 和 
H15915 (5’ -CTC CGA TCT CCG GAT TAC AAG 
AC-3' ) (Xiao et al 2001)， 由 上 海 生 工 生物 技术 服 
务 有 限 公 司 合 成 。PCR 反应 在 Eppendorf 
Mastercycler Gradient PCR 仪 上 进行 ， 反 应 体系 为 
25uL: 10xbuffer 2.5uL, dNTP lnL (2.5mmol/L), 
引物 L14724、H15915 各 1uL (10nmol/L)， 模板 
DNA 2uL ( 约 50 一 200ng)，0.5nL Tag DNA 聚合 酶 
(2UhL， 购 目 天 根 生物 科技 有 限 公 司 )，17hL 纯 
净 水 〈 购 自 华 舜 生 物 工程 有 限 公 司 )， 加 入 10hL 石 
背 油 。 循 环 程序 为 : 95SC 预 变性 4min; 94'C 35s， 
52C 40S，72C 70s，35 个 循环 72C 延 伸 10min。 
取 4uL 扩 增 产物 经 0.8% 琼 脂 糖 凝 胶 〈 含 EB 
0.5ung/mL) 电泳 ， 电 泳 缓冲 液 为 SxTBE (用 时 稀释 
10 倍 )， 电 压 SVcm， 电 汶 1h 后 ， 用 Gene-Genius 
凝 胶 成 像 系 统 拍照 检测 ,目的 基因 用 UNIQ-10 柱 式 
DNA 胶 回 收 试剂 盒 〈 上 海 生 工 技术 服务 有 限 公 司 ) 
纯化 回收 ， 送 上 海 生 工 技术 服务 有 限 公 司 ， 由 ABI 
377 自动 测序 仪 进行 正 反 癌 测 序 。 测 序 所 用 引物 同 
PCR 扩 增 反应 。 

1.5 数据 分 析 

1.5.1 RAPD 分 析 经 电泳 获得 RAPD 图 谱 , 在 同 
一 电泳 迁移 位 置 上 ， 有 RAPD 扩 增 条 带 的 计 为 1， 
没有 的 计 为 0， 每 群体 所 有 个 体 在 全 部 位 点 上 的 谱 
市 式样 (RAPD 表 型 ) 组 成 一 个 (0 或 1) 和 窍 阵 。 用 
POPGENE version 1.3.1 软件 (Yeh et al, 1999) 进 行 分 
析 ， 计 算 多 态 位 点 数 、 多 态 位 点 比例 CP)、Nei 基 
因 多 样 性 指数 〈 刀 )、Shannon 多 样 性 指数 (7) 及 和 群 
体 间 遗传 距离 等 。 用 Arlequin Ver 3.01 软件 进行 分 
子 方 兰 分 机 (analysis of molecular variance, 
AMOVA ) (Excoffier et al, 2005)。 



































1.5.2 Cytb 基 因 序列 的 对 位 排列 “alignment) 使 
用 Clustal X 软件 (Thompson et al，1997)， 并 在 
SEAVIEW 程 序 (Galtier et al 1996) 中 辅 以 手工 校正 。 
采用 Mega 4.0 软件 (Kumar et al 2004) 分 析 核 音 酸 组 
成 及 统计 单 倍 型 、 计 算 转 换 / 题 换 值 CTUMTv) 及 
Kimura 双 参 数 站 传 距 离 (Kimura 2-parameter 
distance， 开 2-P)， 并 用 该 软件 把 核 背 酸 序列 转化 为 
进 基 酸 序 列 ; 同时 采用 Mesga 4.0 软件 的 邻接 法 
(neighbor-joining, NJ) 构建 大 弹 涂 鱼 各 单 倍 型 间 系 
统 发 育 关 系 树 ， 以 至 弹 涂 全 、 大 鱼 弹 涂 介 、 弹 涂 鱼 
为 外 类 群 ， 以 Kimura 双 参数 法 为 稚 代 模型 ， 系 统 树 
分 文 的 置信 度 采 用 目 展 法 (bootstrap analysis, BP) 
重复 检测 1000 次 。 用 Arlequin Ver 3.01 软件 统计 种 
群 核 音 酸 多 样 度 (nucleotide diversity， 如 、 单 倍 型 
多 样 度 (haplotype diversity，h) 及 其 标准 闫 (SD)。 
采用 AMOVA 分 析 方 法 , 以 1000 次 重复 随机 抽样 单 
倍 型 重 排 后 进行 显著 性 检验 ， 用 于 佑 计 大 弹 涂 鱼 种 
群 遗 传 结构 及 不 同 地 理 种 群 遗 传 变 寞 的 分 布 。 用 
DnaSP Ver. 4.10 软件 (Rozas et al, 2006) 对 群体 间 分 
化 指数 《F,)、 基 因 交 流 值 (NM,) 及 Cyt 5 基因 单 倍 
型 铺 配 进行 分 机 。 为 了 解 大 弹 涂 鱼 的 种 群 历史 ， 以 
Arlequin Ver. 3.01 软件 进行 Tajima’s DD 中 性 检验 。 以 
TCS 软 件 (Clement et al，2000) 建 立 网 状 杀 缘 关系 
(network)。 
根据 公式 7=Dky2w 计 算 谱系 间 分 化 时 间 ,/ 是 核 
并 酸 的 变异 率 速率 (1% 一 2.5%) / 百 万 年 〈Irwin et 
al ，1991; Martin et al ，1992)。 用 公式 t=2ut 佑 算 
种 群 扩 张 时 间 T(Roges et al 1992) ， 7=tx 代 时 , 代 
时 由 研究 对 象 的 性 成 熟 年 龄 决定 。 其中, 1 表示 种 群 
扩张 至 今 所 经 历 的 代数 , t 是 错 配 分 布 的 模型 , 可 由 
Arlequin3.01 软件 分 析 获 得 , yx 是 所 测 Cyt 2 序列 总 的 
变 寞 率 , 用 公式 =24K 计 算 (天 是 序列 中 核 背 酸 的 数 


量 )。 
2 .和 党 


2.1 RAPD 分 析 

本 研究 从 100 个 10 碱 基 随 机 引物 中 ， 算 选 出 
多 态 性 丰富 的 引物 30 个 : S5、S11、S12、S13、S27、 
S28、S37、S43、S45、S47、S102、S105、S107、 
S108、 S112、 S114、 S117、 S119、 S122、 S123、 S124、 
S125、S126、S127、S128、S129、S132、S137、 
S139、S140。 对 于 JDS、DH、XP 3 群体 ， 每 群体 
各 扩 增 20 个 体 。 每 个 引物 都 能 扩 增 出 6 一 13 条 清 
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晰 而 重复 性 好 的 条 带 〈 图 1)。30 个 引物 在 JDS、 
DH、XP 群体 分 别 扩 增 出 236、270、274 条 带 ， 共 
301 个 有 效 位 点 。 每 个 引物 在 不 同 群 体内 平均 扩 增 
出 7.86 一 9.13 条 上 市， 其 中 多 态 性 条 融 达 90% 以 上 。 
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图 1 
Fig. 1 
JDS: 1-10; DH: 11-20; XP: 21-30. 
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DNA 扩 增 大 小 在 0.3 一 3.0kb。 

3 个 群体 大 弹 涂 鱼 的 多 态 位 点 比例 (P)、Nei 
基因 多 样 性 指数 〈 瓦 ) 和 Shannon 多 样 性 指数 (7]) 
分 别 在 93.02% 一 -96.35%、0.3890 一 0.4219 和 0.5618 





Il7 18 19 20 21 22 23 24 23 26 27 28 29 30 





引物 S126 在 大 弹 深 包 3 群体 中 扩 增 结 来 
Amplified result of primer S126 in JDS, DH and XP populations of Boleophthalmus pectinirostris 


表 1 大 弹 涂 鱼 群 体 间 的 多 态 位 点 比例 (P)、Nei 基因 多 样 性 (HH) 和 Shannon 多 样 性 指数 (7) 
Tab.1 Polymorphic loci (P), Nei’s gene diversity (H) and Shannon diversity index (71) among populations 


of Boleophthalmus pectinirostris 


多 态 位 点 数 Polymorphic loci 


多 态 位 点 比例 (P，%) Propotion of polymorphic loci (P，%) 


Nei 基因 多 样 性 (H， 平 均值 + 标准 误 ) Nei’s gene diversity (H, meantSE) 
Shannon 多 样 性 指数 〈I， 平 均值 + 标准 误 ) Shannon diversity index(IL, mean+SE) 


一 0.6044。 由 表 1 知 ， 大 弹 涂 鱼 和 群体 的 遗传 多 样 性 
都 很 丰富 ， 其 中 DH 群体 最 蝇 ， 而 JDS、XP 和 群体 之 间 
比较 接近 ， 却 低 于 DH 群体 。 对 3 个 群体 间 的 遗传 
差 寞 进行 AMOVA 分 析 表 明 ， 总 的 遗传 分 化 指数 五, 
值 为 0.03196， 和 群体 间 瓦 , 值 在 0.03004 一 0.03254 
(P<0.05), 存在 显著 的 遗传 分 化 。 在 遗传 莽 由 的 总 
变异 中 ,有 96.80% 的 差异 来 目 各 群体 内 部 , 仅 3.20%4 
的 差异 来 自 群 体 间 〈 表 2)。3 个 群体 内 的 遗传 相似 
度 都 很 高 ， 在 0.9403 一 0.9553， 而 彼此 间 的 遗传 距 
离 却 不 大 ，JDS 与 DH、XP 间 的 遗传 距离 分 别 为 
0.0457、0.0543，XP 与 DH 间 的 遗传 距离 为 0.0616。 
2.2 Cytb 基因 序列 分 析 

2.2.1 序列 变异 本 实验 中 ， 获 得 的 Cyt 2 基因 的 
核 甘酸 序列 长 度 为 1141bp， 最 后 一 个 碱 基 为 T, 在 
线粒体 上 为 不 完全 终止 密码 子 U (Anderson et al 
1981)。 经 MEGA 4.0 软件 分 析 ， 从 42 尾 大 弹 涂 全 
中 共 获 得 33 个 Cyt 5 基因 单 倍 型 , GenBank 序列 号 
为 EU5$03190 一 EU503222。 每 群体 分 别 有 8 一 10 个 
单 倍 型 ， 其 中 H2、H4、HS5 是 共享 单 倍 型 ， 以 H2 
共享 频率 和 出 现 次 数 最 多 (5 尾 个 体 ， 占 测序 样本 
总 数 的 11.9%)， 其 余 的 30 个 为 各 群体 特有 的 单 倍 























九段 沙 JDS 定海 DH 霞浦 XP 
281 290 280 
93.36 96.35 93.02 


0.3890+0.1421 
0.3018+0.1890 


0.4219+0.1086 
0.6044+0.1418 


0.3891+0.1430 
0.5619+0.1897 





型 ( 表 3)。 碱 基 组 成 中 ，A、T、C、G 含量 分 别 为 
26.1%、32.1%、27.7%、14.1%, A 十 T 舍 量 为 58.2%， 
G 十 C 含量 为 41.8%， 表 现 出 明显 的 反 G 侦 傈 。 在 
1141 个 序列 位 点 中 来 发现 插 入 和 缺失 , 其 中 保守 位 
点 有 1083 个 、 变 寞 位 点 58 个 (图 2)。 序 列 的 平均 
转换 与 颠 换 (TITv〉 比值 为 4.0， 密 人 码 子 第 三 位 点 
的 TV 值 为 7.9， 显 示 变 寞 主要 发 生 在 该 位 点 ， 而 
密码 子 第 一 、 第 二 位 点 的 Ti/Tv 值 较 小 ,分别 为 2.2 
和 0.4。4 个 群体 内 Kimura 双 参 数 (K2-P) 遗传 距 
离 为 0.0030 一 0.0058， 和 群体 间 为 0.0040 一 0.0056， 
二 者 比较 接近 ， 只 有 定海 与 其 他 3 个 群体 间 的 遗传 
距离 都 高 于 和 群体 内 《〈 表 4)。 

在 氨基 酸 水 平 上 ， 若 以 H2 的 序列 为 参照 ， 则 
所 有 单 倍 型 中 发 生变 异 的 位 点 仅 15 个 ， 即 第 3 位 
Ser 变 为 lle 或 Cys， 第 20 位 Asp 变 为 Glu、Asn 或 
Cys, 第 21 位 Leu 变 为 His, 第 22 位 Pro 变 为 Gln， 
第 130 位 Gly 变 为 Cys, 第 135 位 Val 变 为 Ala, 第 
166 位 Gly 变 为 Ser, 第 221 位 His 变 为 Pro, 第 237 
位 The 变 为 Asn， 第 304 位 Val 变 为 lle， 第 307 位 
Leu 变 为 Pro， 第 311 位 Gln 变 为 His， 第 312 位 
Arg 变 为 Pro, 第 376 位 Leu 变 为 Val, 第 377 位 Leu 
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变 为 His 〈 未 列 出 氮 基 酸 序 列 比 较 图 )。 一 0.9818、0.0030 一 0.0058, 平 均值 为 0.9814、0.0048， 
2.2.2 ”种 群 超 传 结 构 4 个 大 弹 涂 鱼 群 体 的 单 倍 型 其 中 XP 群体 呈现 出 最 低 的 单 倍 型 多 样 性 却 显现 出 
多 样 性 〈1、 核 苷 酸 多 样 性 (rr) 值 分 别 在 0.9556 ”最 高 的 核 甘 酸 多 样 性 ( 表 3)。 利 用 DnaSP 估算 大 


表 2 大 弹 涂 鱼 群 体 间 遗传 差异 的 分 子 方差 分 析 (AMOVA) 


Tab.2 Analysis of molecular variance(AMOVA) in populations of Boleophthalmus pectinirostris 





变异 组 成 变异 百分比 
人 ce 自由 度 df | . 
变异 来 源 Source of variation Varlance components Percentage of variance 
Cytb RAPD Cytb RAPD Cytb RAPD 
群体 间 Among pop 3 2 0.053 1.84 1.93 3.20 
群体 内 Within pop 29 57 2.706 55.85 98.07 96.80 
总 变异 Total variation 32 59 2.759 57.69 100.00 
J 


1122333334444445556666666666777788899999999900001111 
96664611125892366998990112466679111812501223356624891123 
788025486987508364508682246256226983057338830016846997911360 

H2 AGGCTCCCTAATCGCATTGGATCAAACGTAACCCCCGAGCCTGTTAGAACACCGTGCT 


H24d:. a Sn eb ne es A 下 


图 2 4 个 大 弹 涂 鱼 群 体 的 33 个 Cyt 2 基因 序列 单 倍 型 的 变异 位 点 比较 
Fig.2 Varlable sites of Cytb sequence among four populaitons of Boleophthalmus Pectinirostris 
图 上 数学 由 上 人 至 下 为 碱 基 变异 位 点 在 所 测序 列 中 的 位 置 。 


Numbers at the top of the figure indicate the variable sites In the sequece analysis based on H2. 
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弹 涂 鱼 群 体 间 分 化 指数 F, 值 与 群体 基因 交流 值 
(CN)， 结 果 表 明 ， 大 弹 涂 鱼 4 个 群体 间 的 Fy 值 仅 
为 0.00043 一 0.07814， 总 的 群体 分 化 指数 玉 , 值 为 
0.01932 〈P>0.05); 和 群体 则 的 基因 交流 值 CN) 均 
远大 于 1 (14.50 一 30.79)， 其 中 ，XP 与 JDS、CX、 
DH 3 群体 间 的 基因 交流 值 (14.50 一 21.65〉 较 小 ， 
而 JDS 与 CX、DH 2 个 群体 间 的 基因 交流 值 〈27.00 
一 30.79) 较 大 ( 表 5)。AMOVA 分 析 的 结果 显示 ， 
98.07% 的 变异 来 目 种 群 内 ， 而 种 和 群 则 的 变异 仪 占 
1.93%( 表 2)。 

以 邻接 法 (NJ) 构建 系统 发 育 树 (外 群 为 青 弹 
涂 鱼 、 弹 涂 鱼 、 大 鱼 弹 涂 鱼 的 Cyt 2 基因 序列 ， 
GenBank 序 列 号 分 别 为 : EU503223 一 EU503225) 
显示 ， 大 弹 涂 但 4 个 群体 可 以 分 为 两 个 大 谱系 
(lineage) 和 多 个 文系 ,文系 A 有 28 个 单 倍 型 构成 ， 
支持 率 为 52%; 支 系 B 由 5 个 单 倍 型 构成 ， 支 持 率 
为 78%。 但 这 两 个 谱系 并 非 以 4 个 群体 按 地 理 单 元 
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分 别 聚 类 ， 而 是 每 个 谱系 都 包括 其 他 群体 的 个 体 
(图 3)。 两 谱系 则 的 序列 拳 寞 (Dxy) 为 0.01， 经 
计算 可 得 7 值 为 0.5 一 0.2 百 万 年 (MY)。 此 期 为 地 
质 时 代 上 更 新 世 晚 期 ， 说 明 大 弹 涂 鱼 在 这 期 间 产 生 
了 遗传 结构 上 的 差异 ， 从 而 呈现 出 现今 cyt 5 基因 序 
列 差 异 较 大 的 两 个 谱系 。 
2.2.3 种群 历史 利用 Arlequin3.01 软件 ， 对 4 个 
群体 分 别 进 行 Tajima’s D 中 性 检验 ， 因 各 群体 样本 
数 较 少 (10 一 11)，Tajima's D 值 均 为 负 值 (-1.542 
一 0.702) 且 未 达到 显著 水 平 (P>0.10)， 故 将 所 有 
个 体 作 为 一 个 整体 进行 分 机， 发 现 Tajima’s D 值 
(CD= -2.111) 达到 显 车 水平 (P<0.05)。 说 明 大 弹 涂 
鱼 曾 作为 一 个 整体 经 历 过 种 群 扩张 〈population 
expansion ) 。 
Network 的 结果 显示 《图 4)， 大 弹 涂 鱼 所 有 的 
Cytb 单 倍 型 拥有 一 个 共享 的 中 心 单 倍 型 CH2)， 并 
聚 成 两 个 分 支 , 其 中 采 自 XP 的 Hl、H26 与 采 自 JDS 




















表 3 大 弹 涂 鱼 群体 样本 数 、 单 倍 型 (H)、 单 倍 型 多 样 性 (hh) 和 核 苷 酸 多 样 性 (和 ) 值 
Tab.3 Number of sample size and haplotypes, haplotype diversity (h) and nucleotide diversity (7) 
of the four Boleophthalmus pectinirostris populations on Cyt b gene 








样本 数 单 倍 型 多 样 性 。” 核 甘酸 多 样 性 
采样 地 Populations Sample 单 倍 型 Haplotypes(H) Haplotype Nucleotide 

size diversity(h) diversity(7) 
上 海 九 段 沙 JDS 11 H2 (2) *, H4, H6, H7, H8, H9, H10, Hl11, H12, H13 0.9818+0.0463 ”0.0053+0.0031 
浙江 慈溪 CX | H2, HS5 (2) *, H14, H15, H16, H17, Hl8, H19, H20, H21 0.9818+0.0463 0.0048+0.0028 
福建 优 浦 XP 10 Hl (2) *, H3 (2) *, H4, H22, H23, H24, H25, H26 0.9556+0.0594 ”0.0058+0.0034 
舟山 定海 DH 10 H2 (2), H5, H27, H28, H29, H30, H31, H32，H33 0.9778+0.0540 0.0030+0.0019 
总 计 Total 42 33 0.9814+0.0118 ”0.0048+0.0026 





H〈) * 内 数字 表示 茶 群 体 中 ， 出 现 该 单 僧 型 的 样本 数 ， 互 表示 共 社 单 倍 型 。 


Figure in the brackets shows the number of speclmen in the population, H Indlicates the shared haplotypes. 


表 4 大 弹 涂 鱼 群 体内 与 群体 间 的 K2-P 遗传 距离 


Tab.4 The average K2-P distance among and within the four 


Boleophthalmus pectinirostris populations 


群体 Pop 九段 沙 JDS ”慈溪 CX 霞浦 XP 定海 DH 
九段 沙 JDS 0.0054 

花溪 CX 0.0051 0.0048 

霞浦 XP 0.0056 0.0053 0.0058 

定海 DH 0.0044 0.0040 0.0048 0.0030 


表 5 


大 弹 涂 鱼 群 体 间 的 基因 交流 值 Nw〈 对 角 线 上 方 ) 与 种 群 分 化 指数 Fst( 对 角 线 下 方 ) 


Tab.S Matrix of pairwise Nm (above diagonal) and F'st (below diagonal) among 
the four Boleophthalmus pectinirostris populations 








Nm 
EE 九段 沙 JDS 。 ”慈溪 CX ”霞浦 XP ”定海 DH 
九段 沙 JDS 27.00 21.65 30.79 
慈溪 CX 0.00043 15.41 30.79 
霞浦 XP 0.0056 0.00513 14.50 
定海 DH 0.04181 0.01102 0.07814 
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图 3 基于 Cytb 基因 序列 构建 的 4 个 大 弹 涂 鱼 群 体 邻接 树 


Fig.3 Neighbor-joining tree based on Cyt gene sequence varlation of four Boleophthalmus 


pectinirostris population 








节点 处 的 数值 为 1000 次 bootstrap 检验 的 文 持 率 ， 仅 显示 文 持 率 大 于 50% 的 数据 。 


Numbers at the nodes indicate bootstrap values with 1000 replicates, only > 50% shown. 


的 H9、H12 及 采 自 CX 的 H15 单独 形成 一 个 分 支 ， 
其 余 28 个 单 倍 型 聚 成 另 一 分 文 ， 但 4 个 群体 的 单 
倍 型 并 没有 依 采 样 地 形成 各 目的 分 文 (clade )， 而 
早 现 出 一 个 非典 型 星 状 (star-like〉 的 发 散 图 。 用 
DnaSP4.0 对 Cytb 单 倍 型 进行 错 配 分 析 , 结果 显示 ， 
所 有 的 Cyt 5b 单 倍 型 呈现 为 一 条 虽 不 完全 平滑 但 也 
只 有 一 个 明显 顶峰 的 单 峰 型 曲线 (图 5)。 这 些 同样 
说 明 ， 大 弹 涂 鱼 历史 上 有 过 种 群 扩 张 。 本 研究 中 ， 
T=2.587，K=1141， 代 时 取 1( 因 大 弹 潜 鱼 通 党 在 1 
龄 即 可 达到 性 成 熟 )， 计 算得 到 了 值 为 0.057 一 
0.023MY， 说 明 大 弹 涂 鱼 和 群体 扩张 时 间 约 为 更 新 世 
晚期 。 











3 讨论 


3.1 大 弹 涂 鱼 群体 的 遗传 多 样 性 

遗传 多 样 性 不 仅 是 形成 生物 多 样 性 的 基础 ， 也 
是 物种 进化 潜能 的 保证 。 遗 传 多 样 性 的 下 降 ， 可 能 
导致 物种 对 环境 适应 能 力 的 降低 ， 这 对 生活 在 对 外 
多 变 环 境 中 的 群体 来 说 是 一 个 极 大 的 威胁 
(Vrijienhoek et al 1994; Frankham et al, 2002)。 以往 
关于 大 弹 涂 全 跟 传 多 样 性 的 报道 很 少 (Jin et al 
2004; Kanemorl et al; 2006; Zhang et al, 2000)，Jin et 
al(2004) 通 过 RAPD 的 研究 结果 显示 ， 浙 江宁 海 三 
门 湾 海 区 的 大 弹 深 鱼 群体 的 壮 传 多 样 性 水 平 比 较 
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低 ， 多 态 位 点 比例 (P) 仅 21%， 呼 吁 加 强大 弹 涂 
鱼 种 质 资源 的 保护 。 但 本 研究 通过 RAPD 技术 及 
Cyt b 基因 序列 分 析 却 显示 ， 长 江口 及 其 南部 邻近 
的 丹 山 、 晨 浦 等 地 区 的 大 弹 涂 鱼 天 然 群 体 的 簿 传 多 
样 性 却 很 丰富 。RAPD 结果 显示 ，JDS、DH、XP 
群体 的 多 态 位 点 比例 CP) 在 93.02% 一 96.35% 间 ， 
远 远 局 于 三 门 湾 海 区 大 弹 涂 鱼 群体 的 21%， 也 比 野 
生 的 小 黄鱼 (Larimichthys polyactis) 群体 的 91.03% 
及 半 滑 碧 饮 (Cynoglossus semilaevis Giintuer)〉 的 
76% 一 80% 高 (Han et al 2007)， 显 示 出 很 丰富 的 遗 
传 多 样 性 。 同 样 ， 本 研究 从 42 尾 大 弹 涂 全 样本 中 
获得 了 Cyt2 基因 全 序列 的 单 倍 型 33 个 , 平均 单 倍 
型 多 样 性 (2) 为 0.9184， 与 东海 和 黄海 的 日 本 饰 
( Engraulis japonicus ) (h=0.958) 相 当 (Yu et al， 
2005)， 表 现 出 丰富 的 遗传 多 样 性 。 可 见 , 无 论 从 线 
粒 体 DNA 还 是 从 核 基 因 组 DNA 角度 ，JDS、DH、 























图 4 _ 4 个 大 弹 涂 鱼 群 体 Cyt2 基因 序列 单 倍 型 间 的 网 
状 杀 缘 关 系 图 
Fig.4 Minimum spanning networks showed genetic 
relationship among haplotypes of Cyt b gene for 
the four Boleophthalmus pectinirostris populations 
H1 一 H33 表示 蛙 倍 型 ， 空 心 小 圆 图 表示 缺失 的 中 间 单 倍 型 。 
了 1I-H33 indicated haplotypes and small empty circles Showed the 
haplotypes that didn’t present in the samples but necessary to link all 


observed ones in the study. 
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图 5 4 个 大 弹 涂 鱼 群体 Cytb 基因 序列 单 倍 型 错 配 分 布 图 





Fig.S$ Observed pairwise sequence difference and the 
expected mismatch distribution under the model of 
population expansion on the Cyt b gene haplotypes 
for the four Boleophthalmus pectinirostris 
populations 


CX、XP 野生 大 弹 涂 鱼 群 体 的 遗传 多 样 性 都 很 高 ， 
具有 较 高 的 进化 潜力 ， 这 为 开展 不 同 地 理 种 群 的 种 
杂交 和 选 育 提供 了 种 质保 障 。 
3.2 ”种 群 遗传 结构 及 种 群 历史 
种 群 分 化 指数 (Fy) 常用 来 表示 两 个 种 群 间 的 
遗传 分 化 程度 ， 在 0 到 1 的 范围 内 ，F 值 越 大 ， 两 
种 群 的 分 化 程度 越 高 (Wright et al 1978)。 种 群 基因 
交流 值 (CN,) 用 来 表示 种 群 间 的 基因 交流 程度 ， 寿 
No>4， 则 表明 群体 间 是 一 随机 交配 的 群体 (Hartl et 
al, 1989)。 本 研究 Cyt b 基 因 序 列 分 析 显 示 ， 除了 DH 
与 XP 群体 间 的 Fst 值 〈0.07814) 达到 中 等 分 化 程度 
外 ， 其 他 群体 则 的 Fs, 值 均 很 小 (0.00043 一 0.04181) 
(P>0.05)， 而 和 V 值 却 远 大 于 4 (14.50 一 30.79)。 从 
群体 间 K2-P 遗 传 距离 来 看 〈 表 4)， 除 DH 与 其 它 3 
群体 间 的 距离 都 局 于 群体 内 外 ， 其 余 群 体 间 的 平均 
遗传 距离 与 群体 内 接近 。AMOVA 分 析 表 明 ， 仅 
1.93% 的 遗传 变异 来 目 群 体 间 。 这 些 结果 均 表明 ， 
本 研究 所 分 析 的 4 个 大 弹 涂 鱼 群 体 间 的 遗传 分 化 程 
度 较 低 ， 各 采样 点 群体 间 基 因 交 流 频 繁 。DH、CX、 
JDS 之 间 的 基因 交流 值 很 大 《〈27.00 一 30.79)， 可 能 
与 这 三 个 采样 地 彼此 间 的 地 理 位 置 较 近 有 关 。 但 
RAPD 分 析 结 果 却 显示 JDS、DH、XP 3 群体 间 存 在 
显著 的 遗传 分 化 (P<0.05), 表明 在 群体 遗传 差异 的 
检测 上 ，RAPD 比 Cyt 2 基因 更 为 灵敏 ，Han et al 
(2007) 也 发 现 同 样 的 现象 Cyt 2 基因 是 线粒体 上 的 
一 个 基因 位 点 ， 它 是 重 白 质 编码 基因 ， 进 化 速率 受 
到 限制 ， 而 RAPD 技 术 分 析 的 是 核 基 因 ， 其 遗传 方 
式 与 线粒体 基因 不 同 ， 且 RAPD 能 对 整个 基因 组 的 
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多 个 位 点 的 超 传 变 开 信息 进行 扫描 ， 这 可 能 是 导致 
两 种 分 析 方 法 得 出 不 同 凌 传 分 化 结论 的 原因 。 

以 邻接 法 CNJ) 构建 的 大 弹 涂 鱼 Cyt b 单 倍 型 
系统 发 育 关 系 树 〈 图 3) 及 网 状 杀 缘 关 系 图 〈 图 4) 
均 表 明 ，4 个 群体 可 聚 为 两 个 谱系 ， 但 这 两 个 谐 系 
与 地 理 分 布 不 相关 ， 表 明 大 弹 涂 鱼 正 处 于 详 系 重 排 
(lineage sorting) 阶段， 这 与 Kanemori et al(2006) 
的 研究 结果 相似 : 他 们 通过 线粒体 控制 区 序列 分 
析 ， 发 现 有 日 本 有 了 明海、 九州 八代 湾 的 大 弹 涂 鱼 群 体 
虽 存 在 一 定 程 度 的 遗传 分 化 ， 但 并 没有 投 采 样 点 的 
地 理 单 元 分 开 。 类 似 的 现象 在 钩 牙 黄石 首 鱼 
(Macrodon ancylodon ) 、 要 友 鲜 (BEipzzaa1osa 
fimbriata) 等 海水 鱼 类 也 存在 (Durand et al，2005; 
Santos et al, 2006)。 全 于 引起 两 个 谱系 群体 大 小 到 
于 的 原因 ， 疝 需 更 深入 的 研究 种 群 历史 来 解释 。 

Tajima's D 中 性 检验 利用 以 推测 种 群 曾经 历 的 
历史 事件 。 如 果 Tajima’s D 呈 负 值 ， 且 在 统计 学 上 
达到 显著 标准 ， 则 说 明 序 列 中 含有 比 中 性 进化 模型 
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